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ABSTRACT 



This paper presents observations with intermediate spectral and spatial resolution along the 
major and minor axes of the starburst galaxy NGC 253. The spectral ranges analyzed are in the 
region of the stellar MgI6 (~ 5175 A), the near IR Call triplet (~ 8550 A) absorption features, 
and the region of the Ha emission line. We have compared the shape of the stellar features 
with those of reference stars and determined the linc-of-sight velocity distribution (LOSVD) 
of the stellar component by using a two-dimensional Gaussian decomposition algorithm, and 
show for the first time the rotation curve of the stellar component in NGC 253. Comparing 
the rcccsion velocity curves of the gas and stars, it is seen that the stellar component has a 
decoupled kinematics respect to the gas, displaying a shoUower velocity gradient and larger 
velocity dispersion that the gas in the inner regions. The minor axis kinematics together 
with the kinematics across the central 40" along the major axis, suggest the presence of a 
rotating body with a kinematically misaligned axis respect to the main disk of the galaxy. The 
asymmetries in the LOSVD along the minor axis together with the steep velocity gradient of 
the gaseous component suggest a merger scenario to explain the these kinematically signatures. 
The enclosed mass in the central regions is computed to be (2.4±0.5) x 10^ Mq for a radius 
of r = 0.7" (10 pc). A double gaseous component in the central 6 arcsecs is detected from the 
[SIII] A 9069 A data along the minor axis; this seems to be the signature of a superbubble due 
to a supernova rate of 0.05 yr~^. 



Subject headings: galaxies: individual (NGC 253) — galaxies: kinematics and dynamics — 
galaxies: spiral — galaxies: starburst — dust, extinction 
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1. INTRODUCTION 

The well known starburst galaxy NGC 253 (e.g. Rieke et al. 1980), classified as type SAB(s)c (de 
Vaucouleurs et al. 1991) has been subject of many detailed kinematical studies of its gaseous component. 
These studies include long-slit optical spectroscopy (Ulrich 1978; Munoz-Tunon,Vilchez & Castaneda 1993; 
Arnaboldi et al. 1995,hereinafter ARN95), which shows large deviations from non-circular motions. These 
deviations have been attributed to the outflow of gas due the presence of massive star formation (e.g. 
Heckman et al. 1990), as well as the presence of heavy obscuration in the central regions (Scoville et al. 
1985; Sams et al. 1994). The large amount of extinction seems to be the reason for the steeper rotation 
curve obtained with near-IR emission lines compared to optical lines (Puxley & Brand 1995; Prada et 
al. 1996). The infrared measurements are in agreement with the ^^CO (J=l-0) velocity field derived by 
Canzian et al. (1988). The most complete analysis of the gaseous component in the innermost regions 
has been done using the radio recombination H92a line by Anantharamaiah & Goss (1996) (hereinafter 
AG96). They determined a velocity field in the central 150 pc displaying a distinct gaseous kinematics 
subsystem, cxhibing rotation in a plane perpendicular to the galactic disk, and an inner region with possible 
countcrrotation in the plane of the disk. The authors suggest that this result may be indicative of a 
secondary bar inside the known primary bar, or alternatively, a signature of a merger or an accretion event 
during the history of the galaxy. 

The analysis and comparison of the stellar kinematics with the gaseous motions would be of great 
interest in determining the nature of this peculiar kinematics found in the central regions of the galaxy. 
However, up to now the information on the stellar kinematics of NGC 253 is poor and only the velocity 
dispersion in the nucleus is known (Lutz & Prada 1995; Oliva et al. 1995). The study presented in this 
paper reports the results of spectroscopic observations taken in the region of the MgI6, the Call IR triplet 
stellar features, and the Hq emission line along the major and minor axes of NGC 253. 

The paper is structured as follows: Section 2 presents the observations and data reduction. Section 
3 analyzes the gaseous kinematics and the starbust in the center of the galaxy, and Section 4 determines 
the stellar kinematics and evaluates the enclosed mass in the nuclear region. Section 5 is dedicated to the 
comparison between the stellar and gaseous kinematics. 
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2. OBSERVATIONS AND DATA REDUCTION 

During several campaigns, we obtained long-slit spectra along the major (p. a. 52°) and minor (p. a. 
321°) axes of NGC 253 with the double arm ISIS spectrograph on the 4.2 m William Herschel Telescope 
at the Observatorio del Roque de los Muchachos, La Palma, Spain. We used a 1124x1124 Tek CCD with 
0.36" pixcl^^ in the red and blue arm of ISIS respectively. The spectral dispersion was 0.78 A per pixel 
(with a FWIIM=2 pixels); CuAr and CuNe lamps were used for spectral calibration. The observations are 
summarized in Table 1. In the observations of November 1995, the slit was 0.9" wide, slightly undersampling 
the ^ 0.9" seeing, and was placed along the major axis centered on the apparent optical nucleus. In the 
observation of July and August 1996 the slit was placed at the dynamical center given by AG96 with a slit 
width of 1.3" and seeing of 1.2". The slit was 4' long in all the runs. In each of these runs, standard stars 
with spectral types ranging from G5III to K4III were also observed. These stars were used as templates 
to determine the stellar line of sight velocity distribution (LOSVD) of NGC 253 as it will be explained 
in next sections. Standard data reduction including bias subtraction, flat field correction and wavelength 
calibration was carried out using the IRAfQ package. The sky lines were subtracted using template spectra 
from the ends of the slit. A distance of 3.4 Mpc (Sandage & Tammann 1975) for NGC 253 was adopted, 
yielding a scale of 16.5 pc arcsec""'^. 

3. THE GASEOUS KINEMATICS AND STAR FORMATION ACTIVITY 

Figure 1 shows the gaseous position-velocity curves along the major axis of the galaxy, determined by 
fitting Gaussian profiles to the [Nil] A 6584 A and the [SII] AA 6717 6731 A Ha emission lines. The radial 
velocities are heliocentric and the zero position for the spatial axis comes from the continuum peak at ^ 
6500 A. The systemic velocity of Vhel =230 ± 10 km s^^ is in agreement with previous measurements (e.g. 
ARN95). The results for the three species are in agreement and show that the center of symmetry for the 
gas rotation curve is shifted 11" SW from the continuum peak as was first reported by ARN95. Figure 2 
compares the [SIII] A 9069 A velocity curve (measured in November 1995) with the near-IR II2 A=2.1218 
fim position-velocity curve determined by Prada et al. (1996), and with recent measurements by AG96 

''IRAF is distributed by National Optical Astronomical Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under contract with the National Science 
Foundation 



- 5 - 



using the radio recombination line H92a. There is a good agreement between the three species except for 
the SW of the nucleus suggesting that this side is especially dusty. This velocity difference has been also 
found by Prada et al. (1996) comparing near-IR data with optical data, and was attributed to an optical 
depth effect due to dust extinction (see Section 5). 

Figure 3 compares the minor axis position-velocity curve of the [SIII] emission line with the Hq data 
by Ulrich (1978) and Heckman et al. (1990), as well as for the measurements by AG96 using H92a. The 
[SIII] spectra show two components in the central 6": the brighter component has the systemic velocity 
(Viff;L=230 km s~^) while the fainter one is blueshifted 140 km s~^. Further out from the center, our data 
agree with the Ulrich (1978) data, whilst the data by Heckman et al (1990) given only for the SE, show two 
components beyond 10" from the center. We believe that the deeper penetration through the dust of the 
near-IR wavelengths allow us to see into the equator of the galaxy the two components that Heckman et. 
al (1990) and Schulz & Wegner (1992) attribute to a starburst-driven superwind at the surface of a cone 
or bubble. However, our [SIII] data do not agree with the very steep H92a velocity gradient reported by 
AG96 in the inner 4" (see Figure 3); this will be discussed in detail in Section 5. 

Figure 4 presents the velocity dispersion curve along the major axis of NGC 253 for the optical and the 
[SIII] A 9069 A emission lines. The optical velocity dispersion curve shows a depression with a minimum 
value of 65 km s~^ at the dynamical center (11" from the continuum peak), and reaches ^ 100 km s~^ at 
10" at either side of this. For the [SIII], the velocity dispersion reaches a value of 85 km s~^ at 3.5 SW 
from the continuum peak and decreases slowly towards larger radii. The value of the velocity dispersion at 
the dynamical center is in agreement with the measurements by Oliva et al. (1996) using Br7 for a central 
aperture of 4"x4". 

The two components present in the [SIII] emission lines in the central 4" along the minor axis of 
NGC 253, can be interpreted as a large-scale expansion or superbubble due to supernova (SN) events. The 
measurement of the gas expansion velocity at the position where the slit interset the disk of the galaxy, 
namely the expansion velocity of the superbubble, can be used to constrain the SN rate in the nuclear 
regions. Colina & Percz-Olca (1992) consider a SN rate in the nuclear region of 0.05 yr~^, in agreement 
with the value obtained by Heckman et al. (1990). Taking our measurements, the expansion velocity is ~ 
75 km s~^ deprojected using an inclination of 78.5° (Pence 1981) and the density of the medium in the 
bubble, no=630 cm~^ (Heckman et al. 1990), we estimate a radius of 130 pc for the superbubble in the 
disk of the galaxy. The radius can be estimated from the equations given by McLow & McCray (1988) 
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for the mechanical luminosity of a superbubble Lsn ~ 3.6 x lO^^no {cm~^)R^{pc)V^{km s~^)^. As the 
superbubble should be located interior to the edge-on nuclear ring of 20" (330 pc) in diameter reported by 
ARN95, this introduces a constrain in in the SN rate with values much smaller than the one derived by by 
Rieke et al. (1988) and Forbes et al. (1993) who estimates a SN rate of 0.1 or 0.2 yr~^. For this reason we 
believe that the average SN rate in the central 300 pc is ~0.05 yr~^. The value of 0.1-0.2 yr~^ could be 
more local since it was determined from the presence of few radio sources (Rieke et al. 1988), and the [Fell] 
luminosity in the central 50 pc (Forbes et al. 1993) of the galaxy. Finally, we would like to mention that all 
the authors agree for a SN rate of ~0.3 yr~^ in the nucleus of M82. 

4. THE STELLAR KINEMATICS AND ENCLOSED MASS 

We have determined the LOSVD along the major and minor axis of NGC 253 from the long-slit 
MgI6 (5175 A) and Call IR triplet (8550 A) absorption spectra by means of a two-dimensional unresolved 
Gaussian decomposition algorithm. This method has been applied succesfuUy to determine the stellar 
kinematical structure of NGC 7331 (see Prada et al. 1996). Basically, this approach considers that the 
spectrum of the galaxy in a given angular position is the result of a mean stellar spectrum broadened by 
the internal kinematic structure of the galaxy and by the observational-experimental configuration. The 
intrinsic broadening represents the convolution of this mean stellar spectrum with the LOSVD. Additionally 
the instrumental configuration (slit- width, etc), and the atmosphere (seeing) sliglthly broaden the spectrum 
in the plane of the detector, introducing correlations between neighbouring pixels along both the spectral 
and the spatial axes respectively. We parametrize the LOSVD as a sum of two-dimensional Gaussians in 
the plane velocity-angular position. The dispersion of the Gaussians reflect the resolution in both axes, i.e. 
spectral and spatial resolution respectively. The dispersions used for the observations presented here, were 
^2 pixels and ^1.5 pixels in the spectral and spatial axes respectively, with relative separation 3 and 2 
pixels. The algorithm makes a two-dimensional minimization of the difference between the observed and 
the modelled spectra using a Netwon-Raphson iterative scheme. With respect to the previoTis algorithm by 
Kuijken & Merrifield (1993), the two-dimensional decomposition takes into account the correlation between 
adjacent angular positions mentioned above. The use of this additional information produces little change 
in regions of moderate o high signal to noise ratio (SNR) but might be important in regions far away of the 
galactic center. 

The spectral regions analyzed were 5120 — 5290 A and 8450 — 8700 A along the major and minor axes. 
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As velocity template for the November 1995 observations, the KO giant HR218 star was used. The K4III 
HD 213947 and the G5III HD4388 stars were used for the July and August 1996 runs respectively. Figure 
5 shows the spectra in several positions along the major axis in the region of the stellar Mglb feature. The 
emission feature at ~ 5200 A corresponds to the gaseous line [NI] AA 5198.5 5200.7 A. Also plotted are 
the modelled spectra obtained by convolution of the derived LOSVD with the spectrum of the velocity 
template (excluding the [NI] feature). Several absorption features including the main components of the 
MgI6 are reconstructed indicating the rcability of the algorithm, and that the stellar reference spectra is a 
good representation of the mean stellar spectrum in the galaxy. Figures 6 and 7 show the equivalent plots 
in the region around the Call IR triplet along the major and minor axes respectively. In the case of the 
major axis, the spectra correspond to the August 1996 run. In both cases the Call line at 8662 A has been 
excluded from the fit because it appears strongly contaminated by sky emission. The reconvolved spectra 
are very similar using template stars of types from G5III to K4III. Exhaustive analyses of the performance 
of the algorithm, including Monte Carlo simulations will be presented in a forthcoming paper (Gutierrez & 
Prada 1997). 

The LOSVD along the major axis obtained from the Mglb and Call IR triplet are symmetrical and 
show a single component. Figure 8 shows the radial velocities along the major axis; these have been 
obtained by fitting a single Gaussian to the LOSVD. These radial velocities are heliocentric and the zero 
position have been chosen from the continuum peak in each case. The Mglfe velocity curve goes from —90" 
to +120" from the contimmm peak. Ovcrplotted is a linear fit which reproduces adequately the general 
trend of the curve and which has a velocity gradient of 1.7 km s^^ arcscc^^. Separation of this trend seems 
to occur in the region from ^ +40" to +90". In the case of the Call triplet, the data of August 1996 allow 
to follow the velocity curve in a large spatial range (from —25" to +25" from the continuum peak). Also 
shown is a fit to the Call data of a circular velocity motion (see Bertola et al. 1991). The Call velocity 
curve shows a steeper gradient compared to the Mglb curve for the SW side probably due to the presence 
of dust in the SW side of the nucleus (Scoville et al. 1985; Sams et al. 1994). 

Figure 9 presents the Call IR triplet velocity curve for the minor axis obtained by fitting a single 
Gaussian to the LOSVD. The curve shows a velocity gradient which is not expected from circular motions 
along the main plane of the galaxy. Furthermore, the LOSVD shows some degree of asymmetry towards 
the blueshift in the inner 10", as it is shown in Figure 10. 

Figure 11 shows the stellar velocity dispersion curve from the Call IR data along both axes. The 
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velocity dispersion in the major axis has a value of 90±5 km s~^ at 5" SW from the continuum peak which 
is in agreement with the values that Oliva et al. (1996), and Lutz & Prada (1995) give using the CO 2.3 /xm 
absorption feature for a central aperture of 4"x4". Beyond ~ 10" the velocity dispersion decreases down to 
70 km s~^. The velocity dispersion along the minor axis reaches a value of 110±20 km s~^ in the central 
10". This can be attributed to the presence of two kinematically distinct components in the inner regions 
as revealed by the LOSVD. 

Assuming that the velocity dispersion accounts for all the gravitational support (see McGinn et al. 
1989) and following the notation of Gaffney et al. (1993) the enclosed mass for a given radius is 

N ra^(r)(A + B + C) 
M{r < To) = ^^^-^ ^. (1) 

where A = —dlnn{r)/dln{r) , B = —dlna'^{r)/dlnr, and C = V'^{r)/a'^{r). To determine A we have 
adopted a surface brightness profile fi{r) oc 7--0-8±0-2^ valid for 1" < r < 5" determined from near-IR 
observations (Sams et al. 1994). Hence, we can obtain the flux density F{r) oc r^'^ and therefore 
n(r) oc r~^'^. B and C were determined from the results shown in Figures 8 and 11. Figure 13 shows the 
enclosed mass at different radii: for r = 0.7" (10 pc), we obtained an enclosed mass of (2.4±0.5) x 10^ M0. 
This value is in agreement with the enclosed mass found by Gaffney et al. (1993) for M82: M{r < 7.5 
pc)=(3± 1) X 10^ M0. 



5. THE DICHOTOMY OF THE STELLAR AND GASEOUS KINEMATICS 

The direct comparison of our measurements of the optical stellar MgI6 and gaseous (e.g. Ha) 
position-velocity curves reveals a clear dichotomy between the kinematics of the gas and the stars in the 
central 650 pc of NGC 253 (see Figure 12). The effects of dust extinction in determining the line-of-sight 
velocities are very small for the two spectral lines mentioned above to explain such velocity differences 
between gas and stars (see Prada 1995). However, in the SW side of the nucleus, the heavy obscuration 
(Ay~30 mag) reveals a steeper velocity gradient for the stellar Call near-IR feature, which is shallower 
than the Ha velocity gradient, the [SIII] line, the near-IR emission lines (Puxley & Brand 1995; Prada el 
at. 1996), and the radio emission lines (AG96) in this region, but steeper than the MgI6 velocity curve. 

Adopting the larger wavelength velocity data, i.e. the Call IR triplet absorption feature and the 
emission H92a/H2 A2.13 /xm lines as the intrinsic (i.e. extinction free) stellar and gaseous kinematics of 
NGC 253, we conclude that the motions of the stars in the central regions are not coupled with the gas. 
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We think that the velocity difference could be a consequence of the distribution of the gas in a ring. This 
seems to be confirmed by the high resolution measurements made by ARN95 who claim the existence of 
an edge-on nuclear ring due to the presence of an Inner Linblad Resonance at 300 pc from the center. 
They believe that the bar is the responsible for the gas fueling of the nuclear star burst, and therefore will 
influence the dynamics of the gas in the central regions. This suggests that the gas is distributed in a ring 
with ongoing star formation and that the stars will trace the kinematics of the regions inside the ring. This 
is reinforced by the larger value of the stellar velocity dispersion with respect to the gas: 90 and 65 km 
s^^ respectively. This difference also occurs is other galaxies as reported by Terlevich et al. (1990) from 
observations of the Call IR triplet. They find that the only objects which have emission lines substantially 
narrower than absorption lines are starburst galaxies. Similar values than for NGC 253, have been found 
in the nucleus of M82 (Gaffney et al. 1993); in this galax;y, also the slope of the stellar velocity curve for 
the central regions is shallower than the ionized gas as reported by Gaffney et al. (1995). These authors 
conclude that the stellar peak in M82 may be an older bulge population (Gaffney & Lester 1992; Gaffney 
et al. 1993), which might also be the case of NGC 253 given the remarkably similarity with the dynamics 
of M82. 

The kinematics along the minor axis of NGC 253 reveals the presence of an ordered motion which 

appears well described by a circular motion with AV„iarE^40 km s^^ (see Figure 9). The minor axis 
kinematics together with the kinematics across the central 40" along the major axis which has AVma2:^70 
km s^^ suggest the presence of a central rotating body with a misaligned axis respec to the main rotation 
axis of the galaxy. According to the shape of the LOSVD, this body might contains a decoupled kinematical 
component in its center probably due to the accretion of material, either in the form of stars or gas which 
tiggered the starburst. 

The counterrotation reported by AG96 in the inner 2" along the major axis is also present in our [SIII] 
data, although there is a velocity discrepancy in the SE regions due to dust extinction as discussed in 
Section 3. The fact that our [SIII] data match the H92a velocity data along the major axis, suggests that 
dust extinction is not very high in the inner 4" (see Figure 2). Therefore, we do not understand, if this is 

the case, the very step velocity gradient found by AG96 along the minor axis. Furthermore our [SIII] minor 
axis data show two velocity components. We believe that the discrepancy between the two data sets can be 
interpreted as follows: by fitting a single Gaussian to the [SIII] emission lines we can match the H92q: data 
by shifting 3.5" NE our data. This would suggest that the H92a velocity gradient reported by AG96 might 
represent an artefact due to the outflow of gas. On the other hand, it could well happen that we placed the 
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slit not exactly at the dynamical center, since our positioning error was of 1.25", and the rapid rotation 
seen in the H92a data is a kinematic subsystem that might have the same origin than the structure seen in 
the stellar LOSVD along the minor axis. 

This research is based on observations obtained at the WHT, operated by the ING at the Observatorio 
del Roque de los Muchachos at La Palma of the lAC. We thank E. Perez and R. F. Peletier for helpful 
discussions. The spectroscopic observations were performed on service nights in November 1995 and July 
and August 1996. 

This paper is dedicated to the memory of C. D. McKeith, who contributed much to this work before 
his untimely death on Summer 1996. 
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FIGURE CAPTIONS 

Fig. 1. — The optical gaseous position- velocity curves along the major axis of NGC 253, for the Ha (bottom), 
[Nil] A 6584 A (middle), and the [SII] AA 6717, 6731 A (top) emission lines. In the middle pannel the relative 
intensity of the continuum around 6500 A along the major axis is also plotted. 

Fig. 2. — A comparison among the position-velocity curves along the major axis of NGC 253 of the [SIII] 
A 9069 A, the near-IR H2 A=2.1218 /im determined by Prada ct al. (1996), and the measurements by 
Anantharamaiah & Goss (1996) using the radiorecombination line H92a. The dashed lines represents the 
relative intensity of the continuum around 8750 A. 

Fig. 3. — The position-velocity curve along the minor axis of NGC 253 for our [SIII] A 9069 A data,the 
Ha data by Ulrich (U78) and Heckman et al. (HAX90), and the H92a by Anantharamaiah & Goss (1996). 
The central region is plotted magnified on the right. When present, the filled symbols refer to the fainter 
component. The circles in the magnified plot represent a single Gaussian fit to the [Sill] emission lines. 

Fig. 4. — The velocity dispersion curve along the major axis of NGC 253 for the Ha (top) and the [SIII] A 
9069 A (bottom) emission lines. 

Fig. 5. — The spectra of NGC 253 along the major axis in several positions in a region centred at 5200 A. 
The main absorption line is the Mgl6 features. The emision at ~ 5200 A corresponds to the [NI] doublet 
gaseous line. The heavy lines show the reconvolved spectra obtained using our algorithm (see the text). 

Fig. 6. — The same than in the previous plot but in the region of the Call IR absorption triplet (~ 8500 A). 
Only the two Call features at shorter wavelengths (heavy lines) were used to determine the LOSVD. 

Fig. 7. — The same than in the previous plot but the Call velocities along the minor axis. 

Fig. 8. — Heliocentric radial velocities along the major axis of NGC 253 obtained from the Mglfe (bottom), 
and the Call IR triplet (top) absorption features using the method described in the text. The solid lines 
represent a circiilar motion fit to the call data and a linear fit to the Mgl data respectively. The dashed 
lines represent the relative intensity of the continuum at ~ 5000 and 8750 Arespectively. 

Fig. 9. — The same than the previous figure but for the Call data along the minor axis. 

Fig. 10. — The stellar LOSVD along the minor axis in the central regions of NGC 253. Note a certain degree 
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of assymetry towards the blueshift. 

Fig. 11. The steUar velocity dispersion along the major and minor axis of NGC 253 (top and bottom 
respectively). The dispersions were obtained according the LOSVD determined from the Call absorption 
features (August 1996 run). 

Fig. 12. — A comparison between the velocity-position curves for the gas (thin line) and stars (thick line) 
along the major axis of NGC 253 obtained from the analysis of the Ha and Call lines respectively. The 
dashed line is a linear fit to the Mglfe stellar curve. The zero position in the spatial axis corresponds to the 
dynamical center. 

Fig. 13. — The enclosed mass in NGC 253 as a function of the radius. The points and the error bars represent 
an estimate assuming that the velocity dispersion accounts for all the gravitational support. The solid line 
is a linear fits. 
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Table 1: Observations 



Date 


Exposure time (s) 


P A. (°) 


Spectral range (A) 
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Template stars (Spectral Type) 
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